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Abstract 

The Cauchy problem is revisited for the so-called relativistic Vlasov-Poisson 
system in the attractive case, originally studied by Glassey and Schaeffer in 
1985. It is proved that a unique global classical solution exists whenever the 
positive, integrable initial datum /o is spherically symmetric, compactly sup- 
ported in momentum space, vanishes on characteristics with vanishing angular 
momentum, and its £^ norm is below a critical constant Cp > whenever 
f3 > 3/2. It is also shown that, if the bound Cp on the £^ norm of fo is 
replaced by a bound C > Cp, any f3 G (l,oo), then classical initial data exist 
which lead to a blow-up in finite time. The sharp value of Cp is computed for 

all (3 G (1,3/2], with the results Cp = for (5 G (1,3/2) and C 3/2 = § (jf) 1/3 
(when ||/o||£i = 1), while for all (3 > 3/2 upper and lower bounds on Cp 
are given which coincide as (5 { 3/2. Thus, the £ 3//2 bound is optimal in the 
sense that it cannot be weakened to an £^ bound with ft < 3/2, whatever 
that bound. A new, non-gravitational physical vindication of the model which 
(unlike the gravitational one) is not restricted to weak fields, is also given. 
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1 Introduction 



In |GlSc85j . Glassey and Schaeffer inaugurated a series of studies (see [GlScOl] and 
the references therein; see also |HaRe07| ) of what they sanctioned the "relativistic 
Vlasov-Poisson system" (rVP in the following). The rVP poses a classical Cauchy 
problem for a relative density function f t : M 3 x M 3 — > M + of an iV-body system with 
Cauchy data f e(^n (^(dpdqOQ given in form of the kinetic equation 



(d t + v ■ V q + crV q Ml) ■ V p ) ftip, q) = 0, 



'1} 



in which the velocity v G M 3 and momentum p G M 3 of a (point) particle of unit 
mass are related by Einstein's formula (with the speed of light c = 1), 

P /n\ 



the scalar field 6* : R 3 



I- satisfies the Poisson equation 



A q (p t (q) = 4vr / f t (p,q)dp 
Jr 3 

with asymptotic conditional 



i-i 



(3) 



(4) 



Mq) x -\q\ 

when \q\ — > oo, so thatjH (p t = —\id\~ 1 * J ftdp; and a G { — 1, +1} decides whether the 
gradient force field of is attractive (a = —1) or repulsive (a = +1). In this paper 
we are primarily interested in the attractive case a = — 1. 



1 By *P(dpd(7) we denote the probability measures, by ^P„(dpdq) those having finite n-th moment, 
and by (?p„ R £ Q )(dpdq), a > 1, respectively (^p„ n £ 1 )(dpdq), those of these measures which are 
absolutely continuous w.r.t. Lebesgue measure dpdq on R 3 x M 3 (momentum x physical space) with 
density in £ a (dpdg) , respectively in C 1 (dpdq) which are the functions with one continuous classical 
derivative (abusing notation, we identify these measures with their densities). While the relative 
density function f t thus fulfills the requirements of a probability density function, it should really be 
thought of as a continuum approximation to a merely normalized (i.e. relative) empirical "density" 
on (p, <7)-space of an actual individual iV-body system. Incidentally, the conspicuous absence of 
N in dU) — (01) means that we study the evolution on suitable time and space scales. The scaling 
transformation t \— * N~ 1 t, q i— > N~ 1 q, so that v i— » v and p i— > p, together with / i— > N 3 f and tf> t— * <j> 
restores N explicitly in (TIJ , ^ , (gj . Note that (*p n £ 1 )(dpdg) is invariant under this scaling map. 

2 In principle, asymptotic conditions other than ([4]) can be imposed, for instance other harmonic 
behavior indicating "system-external sources at infinity." 

3 Note that cf>t — — |id | _1 * J ftdp does not represent dynamical degrees of freedom beyond those 
of ft. Thus, we will speak of solutions ft of ([T])-((4]). 
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The rVP system is not truly relativistic in the sense of proper Lorentz or even 
general covariance. Yet for the mathematically special (and physically idealized) 
situation of spherical systems, rVP can actually be obtained from truly relativistic 
(and physically relevant) Vlasov models in certain limiting physical regimes; hence, 
rVP may have some physical significance. Its version with o = +1 (denoted rVP + ) 
is obtained directly from the relativistic Vlasov-Maxwell system (rVM) for a single 
specie of electrically charged physical particles with spherically symmetric initial data 
without any further conditions; cf. |Hor90j . Thus, in the repulsive case 4>t(q) can 
be thought of as Coulomb's electrical potential at the space point q at time t. The 
version with a = — 1 (denoted rVP - ) ought to obtain in a "weak field limit" of the 
physically relevant general covariant Vlasov-Einstein system (VE) with spherically 
symmetric data, though we are only aware of some work (see |Ren94] ) on the com- 
bined weak field plus low velocity limit which leads to the familiar Vlasov-Poisson 
system with a = — 1 (VP - ), formally obtained from ((!])- (jlj) by replacing Einstein's 
formula ([2]) with Newton's v = p. Thus, in the attractive case 4>t(q) may be thought 
of as Newton's gravitational potential at the space point q at time t. Unfortunately, 
this gravitational interpretation has to be taken with a grain of salt, for some math- 
ematically interesting phenomena such as stationary bound states [Bat89l IHaRe07] 
and finite-time blow up [GlSc85j (signaling gravitational collapse to a singularity, 
see |LMR08b] ) of this (both psychologically and "physically") attractive version of 
rVP occur in the strong field regime, i.e. were rVP - can no longer be expected to 
be a legitimate approximation to VE. While this would seem to make mathematical 
studies of the strong field regime of rVP - questionable from a physical perspective, 
in the appendix we give an unconventional (and presumably surprising) physical in- 
terpretation of rVP~ with spherical symmetry in terms of distributional solutions 
of rVM for a neutral two-species plasma which is not restricted to weak fields, and 
which could have interesting applications in space physics; thus the mathematically 
rigorous vindication of this electrical interpretation of rVP~ is an important open 
problem. 

In the main part of the present paper we revisit the questions of global existence 
and uniqueness versus finite time blow-up of solutions to rVP~, which were addressed 
already by Glassey and Schaeffer |GlSc85j . We restrict our discussion to classical 
solutions, but ask for the optimal — in the sense of weakest - - constraints that 
guarantee that classical data will launch a unique global solution of the dynamical 
system. In this vein, we prove the following result about rVP - : 
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Theorem 1.1. A unique global classical solution of rVP~ exists for all spherically 
symmetric initial data fo £ H C 1 (dpdq) which are compactly supported in momen- 
tum space, vanish for p x q = 0, and satisfy H/0II3/2 < C3/2; with C 3 / 2 = § (jf) 1 ^ 3 ~ 
0.367. The £ 3 / 2 bound C3/2 on f is optimal in the sense that initial data fo ex- 
ist which satisfy all the hypotheses except that ||/q|| 3 / 2 > C 3 / 2 , and which lead to a 
blow-up in finite time. 

1 /3 

Remark 1.2. The critical case H/0II3/2 = | (if) n °t covered by our theorem. 

Remark 1.3. All global-in-time solutions covered by Theorem X 1.1\ have positive, those 
that blow up in finite time non-positive energy. Among the data that lead to finite 
time blow-up there are indeed some with zero energy. This improves on Glassey- 
Schaeffer's result that negative energy data will lead to finite time blow-up. 

Remark 1.4. By the interpolation inequality, f £ ^ n £^ with (3 > 3/2 implies 
f £ P 1 n £ 3 / 2 , with H/oHp < Clf- 1 '^ implying \\f Q \\ 

3/2 < ^3/2- This shows that 
a global existence and uniqueness theorem analogous to Theorem \1.1\ can also be 
stated with the sharp £ 3//2 condition on f replaced by a sharp £^ condition on fo for 
any (3 > 3/2, with C3/2 replaced by a corresponding sharp constant Cp > ■ 
Beside the sharp C3/2 given in Theorem \l.l\ and the lower bound on Cp for (3 > 3/2 
just stated, we will also give an explicit upper bound on C@ for (3 > 3/2. 

Remark 1.5. The reverse to the interpolation estimates of course is not true: fo £ 
^Pifl£ 3 / 2 with H/0II3/2 < C3/2 does not imply any bound on \\fo\\p for (3 > 3/2. Thus, 
our £ 3 / 2 condition is weaker than any of the possible £P conditions with (3 > 3/2. 
In fact, our £ 3//2 bound is the weakest possible £^ condition for which an analog of 
Theorem \1.1\ can be formulated, in the sense that our £ 3//2 bound on fo cannot be 
replaced by an 2P bound with (3 £ (1,3/2), whatever that bound. Indeed, among the 
data fo satisfying any such 2? bound with [3 < 3/2, there are some with negative 
energy, which lead to a blow-up in finite time by Glassey-Schaeffer's blow-up theorem 
(evidently, \\fo\\ 3 / 2 > C3/2 f or those data). 

Remark 1.6. By the previous two remarks, everything else being equal our sharp 
£ 3 / 2 condition on f improves on the (nonsharp) £°° condition on f in \GlSc8S$ . 

Remark 1.7. // the normalization ||/o|li = 1 is changed to any other value for the 
£ x norm of fo, the values of the critical Cp for (3 > 3/2 change by a simple scaling 
transformation. 
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The rest of the paper is structured as follows. In the next section, we list the 
familiar conservation laws and the virial identities for rVP~. Then, in section 3, we 
find an £ 3//2 -optimal subset of tyi H £ x (dpdg) on which the energy functional of / is 
bounded below, which bound is 0. Section 4 is devoted to obtaining a-priori bounds 
on the data /q. In section 5 we prove our global existence and uniqueness result of 
classical solutions, all of which have strictly positive energy. In section 6, we prove 
that blow-up in finite time occurs for certain data with non-positive energy. We 
comment on the critical case in section 7, and section 8 lists some interesting open 
problems. Finally, in the appendix we give a vindication of rVP~ in terms of certain 
distributional solutions to two-species neutral rVM. 



In our paper we will make use of (many of) the conservation laws, of the virial 
identity, and another identity, all valid a priori for any sufficiently regular solution of 
rVP~. While (most of) these laws and both identities are proved in [GlSc85j (under 
more restrictive assumptions than stated here), for the convenience of the reader, 
these basic results are collected separately in this section. To simplify the notation, 
we will use the abbreviation J for f R3 , and we write 



for the relative density function in physical space. We shall drop the argument dpdq 
from now on from the symbols for the function spaces. 

We begin with the conservation laws for the Casimir functionals of /. Thus, for 
(the pertinent subset of) / G ^3 D £ x we define the g-Casimir functional of / by 



(/) = gof dpdq , for all g : M+ -> R such that g o f G 2} . (6) 



For g( ■ ) = (id (■))", a > 1, we get the a-th power of the £ Q norm of /; when 
a = 1 this yields just the mass functional (= integral) of /. The choice g(-) = 
—id ( • ) log(id (•)//*) gi ves the entropy of / relative to some arbitrary /* G ^3 fl £ 1 , 



Then, since ([T]) is isomorphic to a continuity equation for / on I 6 , we have 



2 Conservation laws and virial identities 




(5) 





(7) 
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Proposition 2.1. Let t i— > / t G ^3 fl be a classical solution of rVP . Then, 
whenever (fo) exists, also (f t ) does, and 

CW(/ t )=CW(/ ). (8) 

Beside the conservation laws just stated, the familiar quantities energy, momen- 
tum, and angular momentum are conserved. 

Proposition 2.2. Let t i— > / 4 G ^Pi fl C 1 fre a classical solution of rVP~ , then the 
energy of f t is conserved, i.e. £(ft) = £(fo), where 



£(f):= J J VT+Wf(P,Q)d P d q -\ JJJJ ^[^^ dpdqdp'dq', (9) 



and the momentum space contribution to £(f) (denoted £ p (f) ) is the kinetic plus rest 
energy, while the physical space contribution (denoted £ q (f)) is the potential energy 
off- 

Moreover, the momentum of f is conserved, i.e. V(ft) = V(fo), where 

V(f):= JJpf(p,q)dpdq, (10) 

and if also f t G ^2 then so is the angular momentum of f , i.e. J (ft) = J(fo), 
where 

J(f) ■= J J ?x pf(p, q)dpdq . (11) 

Beside the angular momentum functional J(f), also the virial functional 

V(/):= J J q-pf(p,q)dpdq (12) 

plays an important role, but it is not conserved. Its time evolution yields what is 
called the dilation identity for rVP~, which in the physics literature would be part 
of a "dynamical virial theorem." 

Proposition 2.3. Let t 1— > f t G H be a classical solution of rVP~ over the 
interval (0,T). Then 

±V(f t )=£ { f,)-JJ 7J L w f< Mdpdq (13) 
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An immediate and entirely obvious corollary of the dilation identity fll3p . which 
nevertheless deserves to be stated in its own right, is the "stationary virial theorem." 



Corollary 2.4. Let 1 1— > f t = / be a stationary solution of rVP~ . Then 

£{s ^ll7^W SMAq - (14) 

Remark 2.5. If the stationary ft = fo has most or all of its mass supported in a 
cylindrical subset of (p, q) space given by B P (0) x IR 3 with P < I, then we can expand 
y/l + \p\ 2 = 1 + \\p\ 2 + 0(|p| 4 ) in (14\ ), both in its r.h.s. and in £(fo), and obtain 
the familiar stationary virial theorem "2E), in = —E pot " of non-relativistic VP~ , viz. 



//| P | 2 /„( P , 9 )d P d 9 =i//^j>d P M,' 



(15) 



Some results will flow entirely from Corollary 12.41 But for the blow-up results we 
also need another part of the virial theorem, identity (26) in [GlSc85j, viz. 

Proposition 2.6. Let t \— > f t G ^3 fl (C 1 be a classical solution of rVP~ over some 
time interval t G (0,T). Then 

^ // \q\ 2 VT+Wft(p,q)dpdq = 2V(f t )- ff \q\ 2 v ■ V,^(g)/ t (p, g)dpdg . (16) 



3 An Xl^-optimal / domain for lower boundedness of £(/) 

The lower boundedness properties of S (/) play an important role in the proof of 
our theorem. Clearly, £ (/) is unbounded below on fl £ 1 , for we can make S(f) 
as negative as we please along the sequence fu G fl £ 1 , given by fnip^q) = 
(167t 2 /9) _1 -R _3 xbi(o)(p) x XB R (o)(q), by letting R j (here, xs is the characteristic 
function of the set S). So one needs to restrict £(/) to some subset of *Pi fl £ 1 . 
Incidentally, though this is not spelled out in the pertinent references, it follows from 
inequalities (5) and (22) in |GlSc85] that £ (f) is bounded below (by 0, then) when 
/ G ^Pi fl £°° with \\fWoQ sufficiently small, and it follows from the displayed but 
unnumbered inequality in the introduction of |HaRe07] that £ (/) > also when 
/ G C\£P with /3 > 3/2 and ||/|L small enough. Neither of these references reveal 
how small is "small enough" and what happens for larger norms or smaller /3. 
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We here are interested in the £^-optimal domain in / space, in the sense that 
for / G fl £^ we seek the smallest possible (3, and the largest possible £r norm 
of /, such that lower boundedness of 8(f) holds, while 8(f) is unbounded below 
when these conditions on / are not met. If the Laplacian in Poisson's equation (J3j) 
for is replaced by the d'Alembertian, then j3 = 3/2 is the optimal (3 value for 
the corresponding "relativistic Vlasov-d'Alembert" (rVdA) energy functional, and 
the sharp value of the critical £ 3//2 norm of / can by computedo Curiously, the 
energy functional 8(f) of rVP~ is bounded below iff the rVdA energy functional is, 
and by the same bound then, even though 8(f) could have been better behaved, a 
priori speaking, because in the rVdA energy functional / and represent independent 
degrees of freedom. To exhibit these aspects very clearly, instead of elaborating on the 
approach of [HaRe07| . which is based on the Hardy-Littlewood-Sobolev inequality 
and standard interpolation arguments and which would also yield the sharp £ 3 / 2 
bound, we investigate the boundedness of 8(f) in terms of the topologically dual 
approach which can be applied almost verbatim also when the Poisson equation is 
replaced by the inhomogeneous wave equation. 

Thus we introduce the following functional of / and 0, 

8(fA)-= jj (VT+W + <f>(q))f(p,q)dpd q + ^ J |v 9 0| 2 (g)d g , (17) 

irrespective of whether the functions / and are related by the Poisson equation 
(J3J) or not, and of whether satisfies (Tj0) or not. If / and are related by = 
— |id| _1 * f fdp, as for solutions to rVP~, then (jHJ) and (IT?]) can be converted into 
one another by using the Poisson equation (J3J) and an integration by parts; in those 
cases 8(f) = min^ g ^i 8(f, ip), and is the minimizer. 

Proposition 3.1. Let Q = {(f,<f>)\f £ 3/2 , G #o}- Then 



inf {8(f,<f))\(f,<j>) G O, \\f\\ 3/2 < (3/8)(15/16) 1 / 3 } 



The £ 3 / 2 bound on f is sharp in the sense that, for any e > 0, C > 0, we can find 
f G H £ 3/2 with \\f\\ 3/2 = | (i|) 1/3 (1 + e), and G 9)q such that 8 < -C. In 

particular, 8 is unbounded below if the £ 3 / 2 bound on f is replaced by any £^ bound 
for any (3 < 3/2. Moreover, the infimum < T73)) is not a minimum. 



4 Unpublished math. phys. seminar talk at the E. Schrodinger Inst., Vienna, Aug. 02, 2002. 
M.K. takes the opportunity to thank N. Mauser for his kind invitation to present these results. 
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Proof of Proposition V3.1\ Using the estimates |p| < a/1 + |p| 2 < 1 + |p|, we have 
K < £ < 1 + /C, where 

W^) : =^ / |V^| 2 (g)dg + |y (b| + 0(g))/(p,g)dpdg, (19) 

so that to prove the lower boundedness vs. unboundedness of 8 (f, 0) as claimed, it 
basically suffices to work with K(f, 0). Only for the precise value of the infimum do 
we need one extra estimate involving a/1 + \p\ 2 . We also introduce the abbreviation 
h(p,q) '■= \p\ + <P{q) for this auxiliary, "ultra-relativistic" single-particle Hamiltonian, 
and we define £~(f) '■= JJ \p\f(p, q)dpdq. 

To prove boundedness below of /C(/, 0) (hence, of £(/, 0)) on the subset of fl 

for which H/H3/2 < § (y§) > we begin by noting that unboundedness below of /C 
can only occur if 0(g) < for some q; hence, we only need to show that the lesser 
functional /C < /C given by 

WA) = ^J |V g 0| 2 (g)dg- J j h4p,q)f{p,q)dpdq (20) 

is bounded below, where h_ := — min{/i, 0} > is the negative part of h. To this 
effect we now apply Holder's inequality to J h_fdpdq, obtaining the estimate 



1 



S7T 



WA) > i~ / |V>| 2 (g)dg- HMJ/IU-n, (21) 



q V\ - ||'"-||rllJ llr/(r-l) ' 



with t still to be determined. A simple integration with spherical coordinates in p 
space gives 

8tt 

TIL# + r) 

whenever ||0_|| 3+T exists. Now the Sobolev embedding says^] i)J(lR 3 ) — > £ a (IR 3 ) iff 

16 
1 - 



\h-\\l=T^T-^U-\\lZ (22) 



a = 6. By (J22J) , this means r = 3, yielding r/(r-l) = 3/2 and \\h-\\l = (tt/15)||0_ "'' 



Using next the inclusion^! supp 0_ C supp 0, we find the estimate 

\h-\\ 3 <(^) 1/3 n\\l (23) 



5 Note that we work with the homogeneous norm; the standard i} 1 embedding of course gives 
^ 1 (R 3 ) £ Q (R 3 ) for all a e [2,6]. 

6 Note that these fine details are necessary only because we did not restrict </> by Poisson's 
equation |3]) with asymptotic condition (j4]); for if we had, then cf> = —(j)- would ensue. 



9 



Hence, we set r = 3 in (I2~lj) and use (l2"B"j) . next recall the sharp Sobolev inequality 
for the embedding fjg — > £ 6 (see |Tal76| ILie83j ). viz. 



V ^ll2- 3 © 4/3 |l^6> ' «21) 



where 3(7r/2) 4//3 is the largest possible coefficient for the ||0|| 6 term, and obtain 



3 / 7T \ 1/3 / 7T \ 1/3 



Thus we have proved that 

inf [s(fA)\(fA) G O, ll/lls/2 < (3/8)(15/16) 1 / 3 } > 0. (26) 

■ 1 /3 

To see that is the best lower bound for £ on f2 when ||/ 1| 3 / 2 < | (y§) , we work 
with trial densities of the type f^{p,q) := ^(p, || 2, with nonpositive e i}o 

still to be chosen. Clearly, > and H/^l^ = 1, and furthermore/C^/^, 0) = 
K.{f 1 0)- We easily compute that 

J j Hp, q)f(p, q)dpdq = -jj^jf (27) 

and ||/*|| 3/2 = II/1-II3/II/1-H2. With the help of (JH and < 0, we now find \\h-\\l = 

(7r/15)||0||g, and we readily calculate that || (j^ = (27r/15)||0|| 5 , so that the trial 

0(g) will have to approach sufficiently fast as |g| — > 00 in order for ||0|| 5 to exist; 
in particular, is fast enough. For such f^, we therefore have 

/ f Hp, Q)f*(p, q)dpdq = -1 |j|jjf (28) 
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and 



1 A5\ 1/3 II0H* 



13/2 2VW ||0||^' 



(29) 



so that 



JjHP,Q)fHp,Q)dpdq = -(^) 1/3 \\f\\ 3/ M\l (30) 
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By taking strong limits we now see that (15U1) is valid for all / G ^Pi fl£ 3 / 2 and G £ 6 . 
In total we have /C(/^, 0) = £7(0), with 

^) : =ilK0||^-(^) 1/3 ||/1| 3/2 ll«- (31) 



!>7r 1 



Now recall that up to translations, the scaling family of functions 

with k > provides us with all the optimizers satisfying (j4]) of the sharp Sobolev 
inequality (|24|) for the embedding i^Q — >• £ 6 , i.e. "=" holds in ff24|) when = K ; see 
[Tal76l[LTe83] . We set = K in note that /*» G^fl so /*» G <Pi H £ Q for 
all a. An easy computation gives us ||0«|| 6 /||0/t|l5 = Sn 1 / 3 /^/ 3 , independent of k. 
Hence, by (J2HD, ||/ 0K || 3/2 = | (§) 1/3 independent of k. Therefore Q (<j> K ) = for all 

k, and since sup| p | >0 (|p| a/1 + \p\ 2 — |p| 2 ) = 1/2, we now conclude that for any small 
e > we have 



£(f"'M = J (VT+W-\p\)f"(p,q)dpdq 
< \ JJ\p\- 1 f >K (p,q)dpdq 

whenever k > C*/e, where is independent of k. The "inf-part" (ITS]) of our 
Proposition 13.11 is proved. 

Incidentally, with a minor modification, the last line of reasoning also shows that 
"inf ^ min." Thus, for any (/, 0) G satisfying ||/|| 3 / 2 < | (y§) , we have 



£UA) = jj(VT+W-\p\)f(P,<l)dpdq + IC(f,<P) 

> J J {^^ + \p\ 2 -\p\)f{P,q)dpdq >0. (34) 

Next, to prove that the £ 3//2 bound on / is sharp, we now show that for any 
given e > and C > 0, we can find / G ty\ fl £ 3 ^ 2 and <f> G S)q such that 
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13/2 = I (ll) 1 ^ 3 (1 + e ) anc ^ ^ — ~~ As stated at the beginning of the proof 
of our proposition, it suffices to prove unboundedness below for K along these lines. 
For this purpose we continue to work with the optimizers of Sobolev's inequality f l24l) . 
<p K , but now invoke the double scaling family of trial densities / = f%\, k > 0, A > 0, 
defined by 

f&(p,q):=i?\*f*{*p,Kq). (35) 

Note that f^ x G n £°° for all k, A, so G <£i n £ Q for all a. Note furthermore 
by an obvious re-scaling of the integration variable that the pertinent density in q 
space is independent of A, i.e. f f%\{p,q)dp = f ft\{P->q)dp f° r a ^ ^> moreover, it 
equals J f^ K \p,q)dp =: p^ K (g), for f^\, K = f <t>K - The density p <t>K (q) = K 3 p* 1 (/tg), 
where p^ 1 is given by 

P 01 (?) = |^(-0i) 5 (9)- (36) 

Thus, not only are f^ K (p,q) and <f> K (q) related by Poisson's equation d3J) with (j3J), 
but so are f%\(p,q) and K (g) for all A. Next, a simple computation yields 

II/?a||3 /2 = ^||/^|| 3/2 , (37) 

with ||/*|| 3/2 = § (i|) V3 , so setting ||/$J| 3/2 = | (±f) 1/3 (1 + e) in ([37} defines a 
branch of a hyperbola k\ = 1 + e in the first quadrant of Cartesian (k, A)-parameter 
space. Along this hyperbola branch, is a probability density with £ 3//2 norm equal 
to (1 + e)Cs/2, i-e. bigger than the acclaimed critical C3/2. Next, a straightforward 
computation shows that along this hyperbola branch, 

k(4W>^) = HIt^ <0 (38) 

for all k > 0. Hence, for any e > and C > there is a unique «(e, C) = 
(32/3tt)(1 + l/e)C > 0, such that iC(f^ {1+e)/K , (j) K ) < -C whenever « > «(e,C). 

Therefore, ^(/wi +6 w sl 0k) is unbounded below, and this implies unboundedness be- 

' — ' 1/3 

low of /C(/, 0) over the set {(/,0) G J], ||/|| 3/2 < C} when G>\ (§) , hence, the 
same holds for £(f, 0) [< 1 + £(/, 0)]. 
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Finally we prove that /C(/, 0) is unbounded below on the domain {0 G #0, / G 
*Pi n ^ ll/lljg < C/3} for any (3 G (1,3/2) and any Cp > 0. (Notice that for (3 = 1, 
which only allows Wf^ = 1, we already proved unboundedness at the beginning of 
section 3.) We ignore ((3D, as allowed by the hypotheses of Proposition 13.11 Thus, 
since for (3 < 3/2 the interval (2,3//?) is not empty, we pick any •& G (2, 3//3) . We 
now choose a family of trial densities of the type f^{p,q) := /^(p, g)/||/i_||^, with 
some nonpositive G satisfying ([3D; will be specified further below. Also for 
this trial family of fs we have }C(f^,<p) = £(/ , 0). Moreover, using again f|22l . we 



find 



ii?+3//3 

= «(''■ ( 39 ) 



and 



1 



K(J*,4) = —\\V q <P\\;-b(0)!^^, (40) 



|2 11^11^+4 
1 0+3 

where a and 6 are some numerical constants dependent on the displayed arguments. 
Now notice that for the stipulated range of (3 and $ values, we have 5 < i? + 3 < 6 
and 5 < "&(3 + 3 < 6, but $ + 4 > 6. Since for G fjg satisfying ((3D we actually have 
G £ a (IR 3 ) for all a G (3,6], it follows right away that H/^H^ < oo. But among 

these functions there are infinitely many for which ||0||^ 4 = oo. Let 0* be such 
a 0. We then can make ||/**||^ < C/3 for any small by multiplying 0* by some 
small constant c if necessary (note that c0* no longer satisfies ((3D, though), while 
||0*||^+4 = oo implies /C(/^*,0*) = — oo. Therefore, we obtain lower unboundedness 
of £(/, 0) on the domain {0 G / G ^ n £^, ||/||^ < C^} for any 1 < (3 < 3/2 
and any Cp. So £ 3/2 in (JTSD cannot be replaced by any £^ with (3 < 3/2. Q.E.D. 

Remark 3.2. Our proof , a variant of which was announced a while ago (see footnote 
4), was inspired by a related proof in [Alytfflj (which in turn was inspired by one in 
[WZS88]) for the energy functional of the nonrelativistic VP system with a = — 1. For 
VP~ the critical (3 = 9/7, and the critical £ 9 / 7 bound on f simply reads ||/|| g / 7 < oo. 

We are now ready to address the boundedness of £(/). 
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Proposition 3.3. We have 

inf {£(f)\f 6 n £ 3 / 2 , ||/|| 3/2 < (3/8)(15/16) 1 / 3 } = 
and the infimum J^Tty noi a minimum. Moreover, we have 



inf {£(/) 



< c 



— oo 



whenever C > Cp, with 

OT) ^ 3(1 ~ 1//3) 



inf 



In particular, Cp = /or a// /3 G (1, 3/2). Furthermore, for all (3 > 3/2 we /iai>e 



15 
16 



1-1//3 



<c«< 



45 



•57T 



5/2 



V/3 



87r2 VnLi(^ + 2/?)r(/3 + 



(41) 
(42) 

(43) 
(44) 



At (3 = 3/2 i/ie upper and lower bounds in coincide; C 3 / 2 = (3/8) (15/16) 1 / 3 is 
the sharp value, i.e. the bound on the £ 3 / 2 norm of f in is optimal. 




C B 0.2 



Upper and lower bounds on Cg together with the vertical line segment between and 
Replacing the co-ordinate label Cp by the diagram acquires the following meaning: 

On and to the right of the vertical line segment, yet on and below the lower curve, we have 
inf £(/) > 0, while inf £(/) = — oo to the left of the vertical line and above the upper curve. 
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Proof of Proposition HOI As to (14"T1) . since the Poisson equation ([3]) together with 
dlj) selects a subset of f2, it follows that £(/) is bounded below whenever £(f,cf>) is. 
Furthermore, since we proved the infimum versus unboundedness of £ (/, 0) when 
(f> E and / G <£i n£ 3 / 2 with pairs of trial functions (/^", </> K ) respectively (/^, 
(with cf) K given by (|32|) ). either pair of which solves Poisson's equation ([3]) together 
with the same infimum versus unboundedness features hold for E (/) when / G 
<Pi n £ 3/2 ; viz. (14 ip is optimal, indeed, and "inf 7^ min." 

As to (g2D and (j43j), we set £(/) := ££(/) + As with £(/, 0) and £(/, 0), 

unboundedness below of K.(f) implies unboundedness below of £(/). In this vein, 
by double scaling we find K.(f K) \) = A _1 £p (/) + K£ g (f), so for any / we can give 
^C(/k,a) an Y value we like by choosing k and A appropriately, at the expense of 
changing \\f\\p to ||/ k ,a|L = (kA) 3 ^ 1 \\f\\p. In particular, once we found a special 
(kq,Ao) pair for which JC(f KQj \ ) < 0, we can then let JC(f K} x) J, —00 by scaling 
along the branch of a hyperbola k\ = kqXq in the first quadrant of Cartesian (re, A)- 
parameter space, keeping ||/ Ki a||^ fixed along this scaling sequence. For each /, the 
borderline case fC(f K> \) = is obtained by choosing k\ = —£^{f)/£ q {f). Hence, 
(1121) holds whenever C > Cp, with Cp the infimum over / G ^1 H £^ of = 

(-^(f)/£ q (f)) 3(1 - m WfWp, i-e. 61. 

As for (JHJ), since (f^\,4> K ) jointly solves Poisson's equation (J3j with (j3J), our 

proof of unboundedness of £(f,(f>) when | ( /"^ a 1 1 3 / 2 — | (if)^ 3 (1 + e ) applies nearly 

verbatim, with 1 1 /"^^ 1 1 3/2 now re P laced h Y \\ft\Wp and I (lf) 1/3 (= II/* 1 II 3/2 = Ca /^ 
replaced by H/^ 1 !!^ = r.h.s. fj44l) . giving (j4*2]) whenever C > r.h.s. ()44p . implying the 
second inequality in (jHJ). The first inequality in (1441) follows from the interpolation 
inequality and (EH), cf. our Remark 1.4. 

Lastly, to prove inf £(f) = —00 when (3 G (1,3/2), we use £(f) = min^ g ^i £ (f, <fi) 

and the fact that inf^ £(/,</>) = —00 for this (3 interval (for which proof we conve- 
niently ignored (j3J) and (T4|)). Yet, it is desirable to also have a direct proof for £(/). 

We need to distinguish between /3 G (1,6/5) and /3 G [6/5,3/2). First, when 
(3 < 6/5, then by taking product densities f(p,q) = g(p)p(q) we have p G £^ for 
/? < 6/5, and then it is well-known that we can find / with £p(f) < 00 and \\f\\p < 00 
but with £ q (f) = —00. Hence Cp = when (3 G (1,6/5). 

Next, when (3 G [6/5, 3/2), we can recycle the trial densities of the type y{p, q) = 
h-{p,q)/\\h-\\^, introduced in the proof of Proposition 13.11 but now with 0(g) re- 
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placed by ips(q) = — e~^/\q\ s , where S > will be specified below; also > will 
need to be re-specified as well. In particular, tps here is not the self-consistent poten- 
tial <p of rVP~, which </> is related to / by Poisson's equation. Rather, ip$ means a con- 
venient way of generating a suitable family of densities p = p^ s = (— V^)^ + Vll^<5||#+3 
which exhibit a local power law singularity; the exponential factor serves to avoid 
integrability problems at spatial infinity that could occur with these power laws. We 
will show that for (3 G [6/5,3/2) we can find 5 > and > such that f*' G 2P 
but G" £ 3 / 2 , and such that < oo while S q (f) = -oo. 

Indeed, since ||/^|L is given by (1551) with replaced by ips, we have J* 5 G" £ 3 / 2 

yet G £^ for any given (3 G [6/5, 3/2) whenever 5 G [2/(2 + 0), 3/(3 + 0/3)), which 
interval is not empty when > 0. For the just established range of 5 values, given 
any j3 G [6/5, 3/2) and for any > 0, we also have p^ s G PI £ 1 . (All of the above 
holds also with (3 G (1, 3/2).) Next, since p G £ 6/,s implies S q (f) > — oo, we certainly 

want p^ s £ £ 6 / 5 , which means we want ||p^ 5 || 6 / 5 = II ^11 6^+3)75/ II V^l #+3 ^° ^ e 00 • 
This gives us the condition 5 G [5/(6 + 2$), 3/(3 + 0)), which interval is not empty 
for any $ > 0. Taking the intersection of our two 5 intervals we find the common 
condition 5 G [max{2/ (2 + 0), 5/(6 + 20)}, 3/(3 + 0/3)), which interval is not empty 
for (3 G [6/5,3/2) when G (0,3/(5/3 - 6)); note that 5/(6 + 20) < 2/(2 + 0) for 
i? G (0, 2) and 5/(6 + 20) > 2/(2 + 0) for > 2. Lastly, we want that < oo. 

Since = W^Oll^ll 

the necessary and sufficient condition for 

£p{f^ 5 ) < oo is 5 < 3/(4 + 0), which may imply, or be implied by, 5 < 3/(3 + 0/9), 
depending on /3; moreover, 3/ (4 + 0) > max{2/ (2 + 0), 5/ (6 + 20)} iff > 2, but in 
that case max{2/(2 + 0),5/(6 + 20)} = 5/(6 + 20). Hence, the resulting common 
intersection of all our 5 intervals is the interval [5/(6+20), min{3/(3+0/3), 3/ (4+0)}), 
which is not empty for (3 G [6/5, 3/2) when G (2, 3/(5/3 — 6)). In summary, for (3 G 
[6/5,3/2), when G (2,3/(5/3-6)) and 5 G [5/(6 + 20), min{3/(3 + 0/3), 3/(4 + 0)}), 
then we have that f^ s G £^ with (3 G [6/5,3/2) but f^ s G" £ 3/2 ; we also have 
£p(f^ 6 ) < oo; furthermore, the pertinent p^ s G" £ 6 ^ 5 . A brief calculation shows 
for 4>(\q\) = — | ici | 1 * J f^ s (p, q)dp that 0'(|g|) x |g| 1 - <5 ( 3+,? ) near the singularity, so 
4> G" H 1 when S > 3/(6 + 20). Thus, 5 > 5/(6 + 20) implies £ q (f^ 6 ) = -oo, and so 

/ £uffi> s \ \3(l-l/(5) 

C/3= P . J ll/^L^O when /3G [6/5,3/2). (45) 



Indeed, JC(f^) = -oo, which implies = -oo. Q.E.D 
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Definition 3.4. In the following, we will call initial data fo E ^Pi fl£ 3 / 2 subcritical 
if II/0II3/2 < C3/2, critical if \\fo\\ z/2 = C3/2, and supercritical if ||/o|| 3/2 > C 3 / 2 . 
We use i/ie analogous terminology for the solutions launched by such data. 



4 A-priori bounds for subcritical / 

By (154"]) . for all critical and subcritical / we have that £(/) = £(f, <t>[f}) > 0, where 

<f>\f] = -\id\~ 1 * J /dp. (46) 

Unfortunately, this does not seem to lend itself to further estimates on the individual 
energy contributions £ p (f) and £ g (f). Yet, for subcritical / we actually have a better 
result than (154]) . and this does lead to a-priori bounds on £ p (f) and £ q (f). 

Proposition 4.1. For any f G ty\ fl £ 3 / 2 satisfying ||/|| 3 / 2 < C 3 / 2; we have 

£{f) a 8^7) / \^n\ 2 m< (47) 

<^3/2 - H/H3/2 



Proof of Proposition 4-1 Since by hypothesis ||/|| 3 /2 is strictly less than the sharp 
critical value, we can retain a little bit from the Dirichlet integral and see, by in- 
specting the steps of the proof of the lower boundedness of £, hence of £ , that now 
we get the estimate ( flTj) . Q.E.D. 

Proposition 14.11 has the following important spin-off. 
Corollary 4.2. Let / G^n £ 3/2 satisfy \\f\\ 3/2 < § (y§) 1/3 . Then, 

\\V g <f>[f]\\l < 8ttx(/)£(/), (49) 
I /" ^TTW/(P, ?)dpdg < (1 + *(/))£(/) • (50) 
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Proof of Corollary The bound (l4"9j) on || V ? 0|| is just a restatement of (HTI) in 
Proposition 14. 1 1 

Noting that for <fi given by (j4*6]) the energy functional ([9]) can be rewritten as 



£(f) = J J VT+Wf(P,q)dpdq-—\\\7 q( f ) [f}\\l (51) 

we see right away that ( l4*9i) now implies ( 150|) . Q.E.D. 

To state our next corollary we need the following lemma, which does not explicitly 
require / to be subcritical. 

Lemma 4.3. Assume that f G ^Pifl£ a for some a > 1 fwft a = oo allowed). Then 
there exists a C(a) that depends only on a, such that the relative density of particles 
(given in (TJPJ satisfies the bound 

IMI 7 <C(«)II/IIS^(/) 1 " ,? , (52) 
where (we recall that) £ p (f) is the kinetic plus rest energy of f (i.e., l.h.s. $5fy)). and 

7 ' 3a- 2' V ' Aa -3 1 j 

Remark 4.4. iVote i/ia£ 77 m ( f53l) zs a decreasing function of a, taking values in 
[1/4,1], while 7 m ( f53j) is an increasing function of a, taking values in [1,4/3]. So 
the optimal possible control of p is with exponent 7 = 4/3, obtained when a = 00, 
while any a < 00 necessarily entails a weaker control on p. In particular, when a = 1 
nothing new is learned beyond what the definition of p says already. 

Proof of Lemma \4-3\ Inspection of the proof, in [Hor81] and [GlSc85] , of the cor- 
responding ||p|| 4/ /3 bound for the relativistic Vlasov-Maxwell equations when / £ 
^Pi R £°° is assumed shows that their proof generalizes to / G ^Pi fl £ a , all a. Thus, 

p(q) = fdp+ f&P 

J\p\<p J\p\> p 

< ( J \f\ a d P y /a (fpY a ' v^+Wfdp 
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which upon optimizing in P yields 

p{q) < C{a){F{f)Y'{GU)f-\ V ■= a/ (4a - 3). (54) 

Raising both sides to power 7, integrating in q and applying Holder's inequality with 
exponent 5 := 3a — 2 = a/^rj), noting that 7(1 — T])S' = 1, yields 

\\p\\l < C(a)-y J(F(f))^(G(f))^ dq 

< C(ay( I (F(/))^dg) V5 ( J (G(/))( 1 -^'dg' 1/5 ' 



C{ay\\f\\ a f( y^TWfiP^^pdq) 



7(1-7,) 



which gives the desired bound by virtue of the hypotheses of Lemma 14.31 Q.E.D. 
As a spin-off of Corollary 14.21 and Lemma 14.31 we now have 

Corollary 4.5. Let f G ^Pi D £ Q for some a > 3/2. Assume furthermore that 

II/II3/2 < I (ll) 1 ^ 3, Then there exists a C(a) that depends only on a, such that the 
relative density of particles satisfies the bound 

\\pl<C(a)\\f\\l(l + H(f)Y^£(fY-^ (55) 



with 7 and r\ related to a by §5, 

Proof of Corollary \4.5\ Under the hypotheses of Corollary 14. 51 we can apply Corollary 
14. 21 which asserts in (!50j) that S p (f) < (1 + x(/))£(/), which bound gives the desired 
bound on ||p|| 7 by virtue of ([52]). Q.E.D. 

Remark 4.6. Note that in 153]) 7 now takes a value in [6/5,4/3], while rj takes 
a value in [1/4,1/2]. For the smallest possible a = 3/2, / f55]) yields only an £ 6//5 
estimate of p, viz\j 

H/4/5 < C||/||$(l + *(f)) 1/2 £(f) 1/2 , (56) 

where C is some numerical constant independent of f. Of course, even weaker £ 7 ' 
estimates of p hold, since under the hypotheses of Corollary \4-5[ 155]) remains true 
if in ( [HP a is replaced by any a' G [1, a), with 7 t— > 7' and 77 1 — > r]' correspondingly. 



incidentally, for d E (0,2) and 5 G [5/(6 + 20), 2/(2 + ■&)), our /> a cx (\p\ - e~^/\q\ s )t G £ 3/2 
but £ £ 6 / 5 . Yet we also have ) = 00 for all 5 G [5/(6 + 20), 2/(2 + 1?)] when # G (0, 2]. 
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5 Subcritical solutions 



We note that the a-priori bounds on the kinetic and potential energy functional and 
on the £ 7 norm of p for subcritical / depend on / only through functionals which 
are conserved by sufficiently integrable classical solutions of the Vlasov evolution. 
This implies at once uniform bounds w.r.t. time on the corresponding quantities for 
sufficiently integrable classical solutions. 

Corollary 5.1. Let t \— > f t G H C 1 , be a classical solution of rVP~ over some 
time interval t G [0, T), and assume that initially and hence for all t G [0, T) we have 

I \ft 1 1 3/2 < I (ll) 1 ^ 3, Then, uniformly in t, we have 

\\VMl ^ 87rx(/„)£(/o) , (57) 
J j + Q)dpdq < (1 + *(/ ))£(/o) , (58) 

and, for all a G [1, oo], 

\\pt\\ y < C(a)||/ ||S(1 + x(/ )) 1 -^(/o) 1 -' 7 (59) 

with 7 and r\ given by 

Proof of Corollary I5.il Under the hypotheses of Corollary 15.11 which imply the 
conservation of energy, S(f t ) = £(fo), and of the £ a norms, ||/t|| Q = ||/o|| Q (note 
that any function / G ^Pi fl is automatically in all £ a ), the uniform bounds 
(EZD (JSBD follow at once from Corollary g~2J and (JSHD from Corollary Q.E.D. 

Theorem 5.2. Let t i— > f t G ^Pi H (t 1 be a classical solution of rVP~ over some 
interval t G [0, T), launched by spherically symmetric data fo with compact support 
in p-space, vanishing for p x q = 0, and assume that initially and hence for all 

1 /3 

t G [0, T) we have ||/t|| 3 / 2 < | (y§) • Then the momentum support is uniformly 
bounded on [0, T) and hence rVP~ possesses a global classical solution. 

The proof of the above, which mimics closely the argument for global existence 
given in [GlSc85], hinges on the following estimate: 

Lemma 5.3. Lett i— > ft be a classical solution of rVP~ over some interval t G [0,T). 
Then for any pair of exponents 7 < 3 < a there exists a constant C aa such that 

|v ? t | < c a ,4f4 e a \\ Pt \\Yp^t) (60) 



20 



where 

P(*):=sup{|p||(p, g )esupp(/ a ), < s < t} (61) 

and 

e:= l^l/3_ e £:=3(1-1/Qt)0 (62) 

1 — 7/a 

Proof of Lemma \5. 3[ We have that t = — |id | _1 * p t and hence 

J\q-q'\<R \Q. ~ Q I ^|?-?'|>-R 19 ~~ 9 I 

< ci(a)|| / o t || aJ R 1 - 3 /° [ + c 2 (7)|| / o t || 7 i2 1 - 3 /T 

< c 1 (a)||/ t || a P(0 3 - 3/a i2 1 - 3/a + c 2 ( 7 )||p 4 || 7J R 1 - 3/7 , 
which upon optimizing in R gives the desired result. Q.E.D. 

Proof of Theorem By spherical symmetry, ft(p,q) = ft(\p\AQ\iP' q) and thus 
Pt{q) — Pt(\l\) so that we have V q 4>t(q) = M(\q\, t)\q\~ 3 q, where 

M(\q\,t) := 4tt f p t (r)r 2 dr. (63) 



Note that limiqi^oo M(\q\,t) = 1. Let the exponent a > 3 be fixed, and let 7 = |^5§ • 
Thus 9/7 < 7 < 4/3. From Lemma [5.31 and Corollary 15. II we then have that in the 
spherically symmetric case 

\V q Mq)\ < mm{\q\- 2 , CI*(t)} < 4(C- 1 / 2 (P(t))^ 2 + |g|)- 2 (64) 

where the constant C depends only on the initial data fo (specifically on its energy, 
its norm, and its £ a norm). 

The next step, following the Glassey-Schaeffer argument, is to analyze the char- 
acteristics of the Vlasov equation in the spherically symmetric case, to conclude that, 
for any T G [0,T), 

P{T ) < P(0) + y/l + P 2 (0) + C 1/2 P(T p . (65) 
We now compute: 

{ = 3 (l-i)i^ = ^|. (66, 
\ ay 1 — 7/a 3a — 3 

Since £ < 2 for a > 3, this implies a uniform bound on P(Tq) regardless of the size 
of C, and since To was arbitrary in [0, T), the theorem follows. Q.E.D. 
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6 Supercritical solutions 



1 /3 

Recall that supercritical solutions are those for which ||/t|| 3 / 2 > g Notice 
that the energy of a supercritical solution can actually be strictly positive, zero, or 
strictly negative, and arbitrarily negative at that. While not much seems to be known 
about supercritical solutions with positive energy, if the energy of a solution is non- 
positive, i.e. if 8 (ft) < (which implies that the solution is necessarily supercritical 
by Proposition I3.3[) . then one can rule out certain classes of global solutions. 

In particular, when £(/o) < we can rule out that f t is stationary. This is an 
immediate consequence of the stationary virial theorem, Corollary 12. 4[ which asserts 
that 8(f) > for any stationary solution. Since stationary solutions are global, we 
thus have ruled out a whole subclass of global solutions when £(/o) < 0. Note that 
the arguments just presented work without any symmetry assumption on /. 

The dynamical virial theorem can be used to rule out further types of global 
solutions ft with 8(f ) < 0. By Theorem III of Glassey and Schaeffer [GlSc85j, for 
8(fo) < no spherical solution with compact support is global. We now prove the 
following generalization of the blow-up result in [GlSc85j . which states that also some 
data with 8(fo) = lead to finite time blow-up. 

Theorem 6.1. Lett \— ► f t E ^flC 1 be a spherically symmetric supercritical classical 
solution of rVP~ over some interval [0,T), and assume that £(/o) < 0. In case that 
8(fo) = 0, also assume that V(fo) < —1/2. Then T < oo. 

Proof of Theorem \6.1\ For £(/o) < this was proved (for compactly supported fo) 
in [GlSc85j . Inspection of their proof of their Theorem III shows that it applies 
verbatim here, where compactness is replaced by the ^3 assumption. 

In case that £(/o) = we argue as follows. By ([13]) we find that V(ft) < 
V(/o) for all t E [0,T), and since by hypothesis V(/o) < —1/2, we conclude that 
V(f t ) < -1/2 for all t E [0,T); in fact, even V(f t ) < -C < -1/2 for all t E (e,T) 
(Note that C will depend on e and on / , but this is irrelevant for the argument). 
Next, as estimated in [GlSc85j in their proof of Thm. Ill, sphericity implies that 
\ff \q\ 2 v ■ V q (j>t(q) ft(p, q)dpdq\ < 1. Thus, by integrating (fTBl) over t from to T, 
and then using the above estimates, we find that 

J I \q\ 2 Vl+WfT(p,q)dpdq <JJ \q\ 2 ^l + \p\^f (p, q)dpdq+2Ce-(2C -1)T, (67) 

and since 2C — 1 > 0, it follows that the r.h.s. of fl67|) < when T is large enough, 
while the l.h.s. is strictly positive. Therefore, T cannot be too large. Q.E.D. 
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Remark 6.2. We don't know whether the condition V(/o) < — 1/2 zs actually sharp. 
Remark 6.3. Classical solutions of nonrelativistic VP" don't blow up \Pfa9%\ [Sch91^ . 

7 Critical solutions 

Our global existence proof of classical solutions to the Cauchy problem for subcritical 
data relies heavily on the a-priori estimates for subcritical solutions obtained from 
a-priori bounds on the individual energy functionals £ p (f) and £ q (f). When / is 

1 /3 

critical so that ||/t|| 3 / 2 = § ({§) , then our Proposition 13.31 still guarantees that 
inf £{ft) = and further that inf £(f) ^ mm£(f) over the set of critical data, but 
these estimates do not seem to produce a-priori bounds on the p-space and g-space 
energy functionals £ p (f) and £ q (f). So no global existence and uniqueness proof for 
critical classical solutions seems in sight. Curiously, no finite-time blow-up result is 
available either. 

8 Some further interesting open problems 

In addition to the wide open questions about critical solutions, there are a number 
of interesting problems that we would like to draw attention to. 

Some of the assumptions in our global existence and uniqueness Theorem 11.11 
are adapted from [GlSc85j . in particular the compact p-support condition and the 
condition that fo vanish for p x q = 0. We expect that the compact support condition 
can be replaced by a weaker one, like sufficiently rapid decay "at infinity," and that 
fo = for p x q = can be dropped. 

Since all solutions for which global existence and uniqueness has been proved 
are subcritical, one may want to consider the reverse question, whether a unique 
spherically symmetric global classical solution necessarily has to be subcritical. We 
suspect that the answer to this reverse problem is negative, but the problem is open. 
If the answer to this reverse problem is negative, it may be forthcoming more readily 
by studying the special subset of global-in-time solutions furnished by the stationary 
solutions to rVP~, of which there are many examples (see |Bat89t IHaRe07j ). 

Then there is the open question of the sharp values of Cp when f3 > 3/2. 

Another interesting question is that of the weakest possible £ Q norm for classical 
data which controls V0. Our global existence and uniqueness proof for Theorem 
11.11 invokes f G £ a for any a > 3, beside subcriticality. We suspect that a = 3 
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may be the critical a value for this question, indeed, but it would be good to have 
a definitive answer. A variant of this question should become particularly relevant 
when one asks for weaker solutions than classical, as e.g. in [dPLi89j. 

Finally, as pointed out earlier, our non-existence result for stationary solutions 
with S (/) < does not assume spherical symmetry. So one wonders how much of 
our other results generalizes to non-spherical solutions. This may seem primarily 
of mathematical interest, for without the sphericity assumption we are in danger of 
leaving the realm of physical validity of rVP~. Yet, since spherical symmetry is never 
a perfect symmetry of nature, it is important to show that significant qualitative 
results do not sensitively depend on having exact spherical symmetry. 

Note added After our paper was accepted for publication we discovered that 
our open question of the sharp values of Cp for (3 > 3/2 is implicitly answered 
by Proposition 1.1 of [LMR08a] which characterizes a subset of the compact Lane- 
Emden polytropes as optimizers in the variational principle (1.14) of [LMR08a], to 
which our variational principle (T43l) in Proposition 13.31 when restricted to (3 > 3/2, is 
equivalent. The quantitative evaluation of Cp, which requires numerical integration 
of the Lane-Emden polytrope equation, was meanwhile carried out by Brent Young 
and will be reported on elsewhere. 
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APPENDIX 

Derivation of rVP with spherical symmetry from rVM 

We begin by recalling the relativistic Vlasov-Maxwell equations for a two species 
plasma, containing a specie of N + positively, and another one of N~ negatively 
charged particles. For simplicity, it is assumed that all particles carry the same 
magnitude of charge and the same mass, as in an electron-positron plasma, and that 
there is overall an even number of particles so that N~ + N + = 2N . Units are 
chosen such that mass and magnitude of charge both = 1. Then, the particle density 
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functions / t (p, q) G ^Pi PI (C 1 satisfy 

«9 t /^(p, q)+v- V q ff(P, q) ± (^(9) + « x fl t (g)) • V p ff(p, q)=0, (68) 

where velocity v and momentum p are related by Einstein's formula (j2J), while the 
electric field £ t (g) G (£ 2 n ^)(dg) and the magnetic field B t (q) G (£ 2 n £ x )(dg) at 
the space point q G M 3 at time t G R satisfy the evolution equations 

a t Bt( ff ) = -V g x£ t (g), (69) 

«9 t £ t (g) = V 9 x B t (q)-4ir I v(u + f t + - v~ fr)(p, q)dp , (70) 

Jm 3 

supplemented by the constraint equations 

V,-B»(«) = 0, (71) 
V 9 -^(9) = 4tt / K/ f + -^/r)(p,g)dp. (72) 

Here, are the relative numbers of charges in the positive and negative specie; i.e.Jfl 
N~ = u~2N and N + = u + 2N. Global well-posedness of the Cauchy problem for 
small initial data is known, and for large data under the additional assumption that 
no singularities occur near the light cone, see [GlSc88l IGlSt861 IGlSt87l IKlSt02j and 
many further references therein. Note that the constraints propagate when satisfied 
by the initial data E (q) and B (q). 

If we set either u + = and v~ = 1 or v~ = and v + = 1, and ignore the 
pertinent / + , respectively /~ equation from the pair (!68l) . the equations (!68l - (!72j) 
of two-species rVM reduce to one-specie rVM. All spherically symmetric solutions of 
repulsive rVP (i.e. rVP + , viz. a = +1 in (JTJ) — d^J) ) also solve this one-specie rVM. 
Indeed, one-specie rVM simply reduces to rVP + in the special situation of spherical 
symmetry, see [Hor90j . The reasons are: (i) all magnetic effects, in particular all 
electromagnetic waves, vanish in a spherically symmetric solution of rVM, rendering 
the charged particle interactions purely electric; and (ii) any two electric charges of 
the same sign repel each other electrically. 

By contrast, spherically symmetric solutions of attractive rVP (i.e. rVP - , viz. 
a = — 1 in (TO)— ([11) ) do not solve rVM for any choice of v + and u~, nor for any 

8 Note that in (|68|) -(l72 j) we have chosen time and space scaled with 2N rather than N as in rVP. 
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other number of electrically charged particle species, so that any attempt to derive 
rVP~ from rVM for two charged species of electrical particles would seem entirely 
misguided. However, we will now argue (convincingly, we hope) that certain families 
of distributional solutions of overall neutral two-species rVM (without spherical sym- 
metry) converge to spherically symmetric solutions of rVP~. Here is our argument. 

Consider now an overall neutral two-species plasma, so N + = N~ = N, i.e. v + = 
v~ = 1/2. We are interested in the distributional solutions representing the actual 
empirical "densities" of the underlying 2iV-body system, which are sums of atomic 
measures (given below). To be able to work with such singular measures, we first need 
to regularize the equations of rVM by convolution with a smooth positive density 
function g £ with £0(3) symmetry and compact support, satisfying J R3 g £ (q)d 3 q — I, 
which will be removed at the end of our construction]^] The regularized rVM readJn 



d t fif (p, q)+v V q fif (p, q)±[(Qe*Et)(q)+vx(Q e * B t )(q)j ■ V p fif (p, q) = , (73) 

d t B t (q)+V q xE t (q) = 0, (74) 

-d t E t (q) + V q xB t (q) = 2vr f v{ 6e * - fr))(p, q)dp , (75) 

V q -B t {q) = 0, (76) 

V q -E t (q) = 2tt / ( Q£ *(fi+-fi-))(p,q)dp. (77) 

Here we introduced the notation [if which could either mean density functions f t G 
^Pi PI H 1 as before, or true measures /if G tyi. In particular, we can allow pi to be 
the empirical relative "densities" on M 6 at time t — of the underlying 2iV-body 
system, defined as follows. Letting particles with even index be positively charged 
and those with odd index negatively, we can identify any point (j)2k, Q2k)k=i Kemn 
(= iV-positive-charges subspace of 2A^-body phase space) with a singular empirical 
relative "density" on R 6 , 

1 - 

A + (p, q) = - S (P ~ P^Mq ~ qsk\ (78) 



N 
fe=i 



9 Note that this regularization breaks the Lorentz covariance; however, we do not invoke any 
Lorentz transformations in our reasoning, and for continuum solutions, when we let g e — > 5, we 
formally recover the Lorentz covariant relativistic Vlasov-Maxwell equations. 

10 The unfamiliar factors 2ir in (IT5|) and ([77)) are a consequence of v~ = 1/2 = v + , which have 
been factored out from under the integrals and multiplied into Att. 
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and each point (p2fc-i, Q2k-i)k=i e ^-odd with a singular empirical relative "density" 
1 r 

A ~fe ?) = ^ S & ~ P^-i)5{q - q 2 k-i) ■ (79) 
fc=i 

We write (pk(Q),qk(ty)1=i e K 12JV if the point in 2iV-body phase space is the initial 
phase point (at time t — 0) of the phase space trajectory 1 i— ► (p k (t), Qk(t))k=i °f our 
plasma, and the corresponding empirical densities (dropping the quotes from now on) 
are denoted Aq (p, g), respectively Af(p,q). Let the //q be given by some Aq (p, g). 
Then these initial empirical densities together with compatible initial data for the 
fields, Eo(q) and B (q), launch a subsequent evolution under (f73]) - (f7T|) for which the 
fj,f are also given by empirical densities, viz. /if(p, g) = Af(p, g), characterized as 
follows. 

Let the evolution equations for the dynamical variables of each particle, i.e. po- 
sition qk{t) and linear momentum Pk(t) be given by 



dt 

dpfc 
dt 



Pk{t) {8Q) 



9k=Qk(t) y/1 + \p k (t)\ 2 ' 

= e k [(g e * E t )(q k (t)) + q k (t) x (g 6 * B t ){q k {t))\ , (81) 

Pk=Pk(t) 

with e k = — 1 if k is odd, and = +1 if k is even. The above equations are the 
Einstein-Newton equations of motion, equipped with the Abraham-Lorentz expres- 
sions for the volume-averaged Lorentz force that acts on each particle. The evolution 
equations for the electric field E t (q) and the magnetic field B t (q) now read 



d t B t (q) + V q xE t (q) = 0, (82) 

- 27V 

-d t E t (q) + V q xB t (q) = 2tt— e k g e {q - q k {t)) q k (t) , (83) 



N 

k=l 



and the constraint equations are 

V q -B(q,t) = 0, (84) 

^ 2N 

E(q,t) = 27r-J2 e kQe(q-qk(t)), (85) 



k=l 
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altogether known as the classical Maxwell-Lorentz field equations; note that 
and fl85|) are just (1751) and f!77|) with /if given by A t . It was proved recently in 
[Ba DiiOH |KoSpOO| that the dynamical equations fl80l) - fl85|) are globally well posed 



as Cauchy problem in convenient Hilbert spaces; see also |KuSpOO Spo04| . Let 



t i— > (pk{t), Qk{t))k=i e K 12JV be the particle phase space trajectory of a global finite 
energy solution to this Abraham-Lorentj 11 ! model fl5U]) - fl5o^ . Then the corresponding 
dynamical empirical densities Af (p, q) satisfy regularized rVM (173|) - (!77|) in the sense 
of distributions. Thus the Abraham-Lorentz model is entirely equivalent to the 
regularized rVM restricted to empirical densities Af. 

We now note a very important point about the relationship of distributional 
solutions of (regularized) rVM and its continuum solutions. Recall that in the intro- 
duction we pointed out that f t "should really be thought of as a continuum approx- 
imation to a merely normalized (i.e. relative) empirical "density" on (p, g)-space of 
an actual individual A^-body system." In this sense, assume that A" is sufficiently 
large so that for both species Aq / 6 n £ : closely in measure (i.e. w.r.t. 
some Kantorovich- Rubinstein distance). Then on a suitably short time scale the 
evolutions of the Af under f[751) - fT77l) will be reasonably closely approximated by 
solutions of (!73l) - (177j) with the initial data Aq replaced by /q GjPi fl XI 1 , and with 
the initial data for the fields, E and -Bo, replaced accordingly] 12 ! On longer time 
scales, various deviations of the rVM evolution for regular initial data will become 
visible, and as we shall see now, under favorable conditions, one of those long-time 
evolutions is captured precisely by regularized rVP~. For continuum solutions we 
can subsequently let g £ — > 6, recovering rVP~. 

To see this, contemplate that both Ag « f$ = f , with f spherical. Then on 
the conventional short Vlasov time scale we will just find non-interacting perfect gas 
dynamics, for — ff =0, then, in fl70|) and (1721) with v~~ = 1/2 = u + . However, if 
all the particles are i.i.d. by fo, then on a longer time scale we should obtain rVP~ 
for ft, and here is why. First of all, by (approximate) spherical symmetry of the 
2A r -body plus field system we expect that magnetic effects can again be neglected, 



11 These semi-relativistic equations of Abraham-Lorentz electrodynamics are a non-Lorentz co- 
variant regularization of formally Lorentz covariant formal Lorentz electrodynamics with point 
charges (which has the unpleasant feature that its formal equations are mathematically ill-defined 
without regularization). The semi-relativistic Abraham-Lorentz model is chosen purely for the 
ease of the discussion. A fully Lorentz covariant regularized Lorentz model is available (see 
ApKiOl, ApKi02. Spo04|) but is considerably more complicated. 

12 For a rigorous proof of this for a scalar caricature of rVM, see [EKR09 . 
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so we set B = 0. But then (1551) . (152"]) . and (1551) are solved by Coulomb's formula 

= -^E e ^ d id i _1 * &( • - *(*))) (?) ( 86 ) 

fc=l 

Evaluating r.h.s. (15T|) with this formula for £7, setting S = 0, we find for particle £, 

N 

e e (g e * E t )(q e (t)) = - ^ ]T {g £ * V, (|id I" 1 * g E { ■ - q 2k (t)))) (©(*)) 

k=l 

N 

+ ^ E (& * V , id I"' * &( • - (<&(*)) (87) 

fc=l 

where we made use of our convention that positive particles carry an even, negative 
an odd index. Now observe that since by hypothesis all particles are i.i.d. by /o, all 
but one of the force terms in (1571) are i.i.d. random variables, the exception being 

(g e * V q (|id I" 1 * g E ( ■ - q t (t)))) (q £ (t)) = 0, (88) 

which states that the Coulomb self-force on the particle vanishes. But then, if I is 
even, a term from the first sum is missing, and if i is odd, a term from the second 
sum is missing. For all the other terms, since all particles are i.i.d., we can use that 
both Af « f t , and paying attention to the correct normalization!^! we now find that 

e e (g £ * E t )(q e (t)) ~ * V 9 (|id T 1 * (g £ * p t ))) (q e (t)) (89) 

which is independent of ef this means that each particle is acted on by a net at- 
tractive central force, for from each particle's perspective the rest of the system is 
singly oppositely charged, as the rest of the system always contains one more of the 
oppositely charged than the equally charged particles; the force is (approximately) 
central by (approximate) spherical symmetry. We therefore introduce 

0t := ~2^& * (lid (90) 

13 Mathematically this is a nice instance where one might be tempted to set two extremely huge 
numbers N and N — 1 equal, but here their difference matters, and N — (N — l) = 1^0 = iV — N. 
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replace eiE t (qi(t)) by -Vi[(gi(t)) for both positive and negative charges in the 
regularized rVM (noting that V„ and g £ * commute) and also replace A t by f t , upon 
which both equations (1731 reduce to the same regularized Vlasov equation 

(d t + v ■ V, - V^(g) • V p ) f t (p, q) = 0, (91) 

in which the velocity v G M 3 and momentum p G M 3 of a (point) particle of unit 
mass are again related by Einstein's formula ([2]), and where (ftl is given by the r.h.s. 
of fl90|) . Note that fl9~Tl) is decoupled from the Maxwell-Lorentz field equations. If f t 
is sufficiently regular, and uniformly so for all e, we may now let q s —>■ 5 and find 
that the resulting f t is a solution to the Vlasov equation 

(d t +v-V q - V q Mq) ■ V p ) ftip, q) = 0, (92) 
where now 

Mq) = -7^(\id\- l *Pt)(q)- (93) 
Clearly, the scalar field <pt satisfies the Poisson equation 

A q( j> t (q) = 4tt-^- / f t (p,q)dp (94) 

with asymptotic condition 

<jH(q) x -(27V|g|)- 1 (95) 

when \q\ — > oo. By a final rescaling of space and time variables we can get rid of the 
factor 1/2N and thus have obtained rVP~. 

Remark 8.1. The i.i.d. assumption on all particles is very important. Indeed, even 
if we merely assume that the particles of each specie are i.i.d. w.r.t. f Q separately, 
then we can still have that particles of opposite species are strictly correlated, viz. 
wherever a positive particle is located, a negative one is, too. This is the classical 
analog of a neutral gas of "positronium atoms," and no electric rVP~~ will result. 
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